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Abstract 

Equilibrium charge state distributions of stable ^^Ni, ^^Co, and ^^Cu beams passing through a lyum thick Mo foil were measured 
at beam energies of 1.84 MeV/u, 2.09 MeV/u, and 2.11 MeV/u respectively. A 1-D position sensitive Parallel Grid Avalanche 
Counter detector (PGAC) was used at the exit of a spectrograph magnet, enabling us to measure the intensity of several charge 
states simultaneously. The number of charge states measured for each beam constituted more than 99% of the total equilibrium 
charge state distribution for that element. Currently, little experimental data exists for equilibrium charge state distributions for 
heavy ions with 19< Z^, <54 (Z^ and Z^, are the projectile’s and target’s atomic numbers respectively). Hence the success of the 

semi-empirical models in predicting typical characteristics of equilibrium CSDs (mean charge states and distribution widths), has 
not been thoroughly tested at the energy region of interest. A number of semi-empirical models from the literature were evaluated in 
this study, regarding their ability to reproduce the characteristics of the measured charge state distributions. The evaluated models 
were selected from the literature based on whether they are suitable for the given range of atomic numbers and on their frequent 
use by the nuclear physics community. Finally, an attempt was made to combine model predictions for the mean charge state, 
the distribution width and the distribution shape, to come up with a more reliable model. We discuss this new ’’combinatorial” 
prescription and compare its results with our experimental data and with calculations using the other semi-empirical models studied 
in this work. 
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1. Introduction 

The charge state distributions (CSDs) of heavy ions often 
have to be considered in accelerator design and in the devel¬ 
opment of stable and radioactive isotope spectrometers. An im¬ 
portant need of information on charge state distributions arises 
with the use of recoil mass spectrometers (such as the electro¬ 
magnetic mass analyser, EMMA [Hi at TRIUMF, Canada, and 
the fragment mass analyser, FMA O at Argonne National Lab¬ 
oratory) to study nuclear reactions in inverse kinematics. These 
systems use electromagnetic fields to filter out the reaction re¬ 
coils from the unreacted beam and rely critically on the fact that 
ions with different m/Q ratios (where m is the mass of the ion 
and Q its charge) follow in principle different trajectories inside 
a field. In experiments with recoil spectrometers the usage of 
gas-cell targets is common. As a beam of recoil products passes 
through the windows of a gas cell it interacts with electrons 
in the material via atomic charge-exchange processes (electron 
captures and losses). A CSD associated with the material and 
thickness of the window is observed for the particles exiting the 
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target. The different charge states of the recoils follow differ¬ 
ent trajectories inside systems employing electric or magnetic 
fields. Due to the typically limited acceptance of systems such 
as magnetic dipoles, only a small number of charge-states (in 
some cases only one) can be transmitted through them and hit 
the detector for a given field setting. Hence, to determine re¬ 
action yields using a recoil spectrometer, one needs to know 
their charge state distribution. Molybdenum (Mo) foils are ex¬ 
tensively used as windows in gas-cells since they offer various 
advantages. The high resistance of Mo under mechanical stress 
(large Young’s modulus lO) allows the achievement of high gas 
pressures using thin windows. Furthermore, Mo may reduce 
significantly the background from fusion evaporation reactions 
compared to a lower-Z foil. 

In the present study we measured the CSDs of stable ^^Ni, 
^^Co, and ^^Cu beams while passing through lyum Mo foils 
(Z^=42). The beam energies were 1.84 MeV/u, 2.09 MeV/u, 
and 2.11 MeV/u for the Ni, Co, and Cu respectively. The re¬ 
sults of this study were used to check the agreement of semi- 
empirical models developed for heavy ions to the experimen¬ 
tal CSDs. In the following we present a brief description of 
the methods for the calculation of CSDs (subsection |1.1| ), the 
formalism for the semi-empirical models considered (subsec- 
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tion |1.2| ), the experimental procedure followed to extract CSDs 
at the University of Notre Dame (Section |^, our results in a 
detailed comparison and analysis of the successes and short¬ 
comings of various semi-empirical models (Section |^, and a 
conclusion and outlook (Section]^. 


1.1. Methods for the calculation of charge state distributions 


The calculation of CSDs for systems with Z>2 is challeng¬ 
ing since, for every ion/target combination of interest, a com¬ 
plete set of electron capture and loss cross sections must be 
known |[4l . Currently, a number of techniques for the calcula¬ 
tion of cross sections (e.g continuum distorted wave approxi¬ 
mation, plane-wave Bohr approximation) and the treatment of 
the electron exchange processes (e.g quasiground-state model, 
three charge-state model, etc.) are adopted by computer codes 
in order to numerically solve the problem. For example, the 
codes CHARGE and GLOBAL have been designed for the 
calculation of equilibrium and non-equilibrium CSDs, at ener¬ 
gies above 100 MeV/u, for heavy projectiles with atomic num¬ 
ber {Zp > 30) in solid targets 0. For intermediate energies 
above 10 MeV/u (and up to about 30 MeV/u), the program 
ETACHA il calculates the evolution of charge state distribu¬ 
tions for ions with up to 28 electrons; the program has also been 
checked at energies around 2 MeV/u ( 71 . Finally, the code CS- 
Dsim has been developed at calculation of equilibrium and non¬ 
equilibrium CSDs for beams and recoils in gas targets |0 . Gen¬ 
erally speaking, for all the codes mentioned above the ability of 
CSD prediction is limited; theoretical calculations on electron 
exchange cross sections don’t have the required accuracy, so 
experimental values must be used instead. Calculated CSDs 
in regions (energies and atomic numbers) where cross section 
measurements have been done or where the adopted approxi¬ 
mations are correct (e.g at high energies) are more likely to be 
accurate than outside of these regions. 

Apart from computational methods such as those described 
above, semi-empirical models are widely used to calculate equi¬ 
librium CSDs for practical applications (e.g on LISE-r-r (3). 
Important parameters of the distributions, such as the mean 
charge state and the distribution width, are given by empiri¬ 
cal formulas based on experimental data. Following that, sym¬ 
metric (Gaussian) or asymmetric functions are used to model 
the distribution and to calculate the equilibrium fractions. The 
main advantage of semi-empirical models is the simplicity and 
speed of the calculations for any system. However, the calcu¬ 
lations are strictly limited to equilibrium CSDs and for systems 
(ion/target) for which experimental data exist. In the energy 
region up to 10 MeV/u for projectile and target atomic num¬ 
ber combinations 19< Zp,Z^ <54 the semi-empirical models 
are expected to more accurately describe the charge state dis¬ 
tributions than the detailed codes described earlier (CHARGE, 
GLOBAL, ETACHA). However, their agreement with experi¬ 
ment has not been tested extensively since very few experimen¬ 
tal data exist for that region. 


1.2. Semi-empirical models for equilibrium charge state distri¬ 
butions: mean charge state and distribution width 

So far there is no quantitative theory to calculate from first 
principles the mean charge state of heavy ions passing through 
solids. However, in the last decades, experimental data on 
equilibrium CSDs have been accumulated in databases, such 
as (TOl . informing the development of semi-empirical calcula¬ 
tions. Difficulties are also encountered regarding the distribu¬ 
tion widths. For high-Z ions, the electron capture and loss cross 
sections change significantly across shell closures giving rise to 
obvious asymmetries in the CSD data (shell effects) ifTTl [TTl . 
For such CSDs the notion of a single distribution width d is 
not valid any more. However, since an analytical calculation 
of CSD for a heavy ion is very challenging, single widths are 
still in use in order to approximately model charge distributions. 
It is common practice to calculate experimental widths assum¬ 
ing a Gaussian distribution (see Eq. below), even for those 
cases for which non Gaussian distributions are expected. Nev¬ 
ertheless, based on the existing experimental data, a number of 
such empirical models have been developed for the estimation 
of equilibrium distribution widths in gaseous and solid targets. 

A number of semi-empirical models from the literature were 
considered in this study, based on their validity in the relevant 
atomic number, energy range, and on their frequent use by the 
nuclear physics community to predict characteristics of CSD. 
In this section we will briefly present the basic formulation of 
each of these models and comment on their performance. 

The most recent semi-empirical model has been developed 
by G. Schiwietz et al. Ga. A many-parameter formula has 
been fitted on data for about 840 experimental CSDs at vari¬ 
ous energies and atomic number ranges for the projectile and 
target. More specifically the data includes CSDs from solid 
targets with 4<Z^ <83 (although more than 40% of the experi¬ 
mental distributions corresponds to C foils), and projectile ions 
with l<Zp <92, at energies up to ~50 Mev/u. The formula 
extracted by the fitting procedure is: 


Zp(8.29x + x^) 

^ “ 0.06/x + 4 + 7 Ax + 

where: 

(Zf-Zp)2 

Cl = 1 -0.26e-^'^“c- 

C 2 = 1 + 0.03vln(Z,), 


( 1 ) 

( 2 ) 

(3) 

(4) 


with the so-called scaled velocity v, given in terms of the pro¬ 
jectile’s velocity Vp and Bohr velocity vb, by: 

V = (5) 

According to their study, deviations from experiment on the 
mean charge-state values, q, are of the order of 2%. For the 
distribution widths, G. Schiwietz et al. in ca proposes a rela¬ 
tion of the form: 

J rry-0 27 ry0.035-0.0009Zp r/ry —XT - 1 / 

d = w[Zp'^-^'Z, " f(q)f(Zp - q)] (6) 
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where w is a scaled width and 


(7) 


for ion species in the range 18< <92 and energies up to 10.6 

MeV/u. According to their improved model: 


fix) = ^ix + 0317P/)lx. 

In (HI, the width w is plotted versus the number of bound elec¬ 
trons, Nb, of the projectile where - q^xp.- According to 

the scaled solid-state data, w^O.l, which is satisfactorily accu¬ 
rate in our present data (Nt ~10). 

In a different approach, (but also based on available experi¬ 
mental data) J.A Winger et al. ca developed a phenomeno¬ 
logical parameterised formula for the mean charge state which 
has the form: 
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q = Zp[l - expiY^aiXf] (8) 

i=0 

where the reduced velocity X is given in terms of the beam’s 
kinematic factor p, as: 

z = piomizPp*^ (9) 

The parameters ai are defined analytically in the way pre¬ 
sented in ca. No additional information is given about the 
model’s performance and agreement with experiment. In the 
same study, a phenomenological parametrization for the widths 
is also derived. According to this model: 


q = Zp[\- Cex/7(-83.275yS/Z°'^’’] 
x[l - expi-l2.905 + 0.2124Zp - 0.00122Z2)]. (13) 

where C=1 for energies >1.3 MeV/u and C=0.9-r0.0769Ep, 
for E^ < 1.3 MeV/u. Their empirical formula for the widths has 
the form: 

d = ^^(0.07535 + 0.19F - 0.2654F2) (14) 

where Y^qjZp. The width model was designed to be more ac¬ 
curate for Zp >54 at energies above 1.3 MeV/u. 

In addition to the above works for the mean charge state and 
width of the charge state distribution, in the present study we 
also considered four more empirical models that provided for¬ 
mulas either for the mean charge state only, or for the distri¬ 
bution width. In particular these are: the work of H. D. Benz 
ll4l for the distribution width, and the works of K. Shima et al. 
oa, To and Drouin d, and A. Leon et al. |[20]| for the mean 
charge state. These models are presented in detail below; H. 
D. Benz, based on Nikolaev-Dmitriev’s work GSl, provides the 
following simple relation for the distribution width d: 

d = 0.21 Zl'^ (15) 


i=2 i=2 

d = expCZ aiilnXy)[\ - expC^Pi^lnZ,)')] (10) 

1=0 1=0 

The parameters ai and Pi are different in this formula from those 
in the previous equation (for details see ref. In order to 

avoid large deviations due to shell effects they used experimen¬ 
tal data at energies above ~5 MeV/u. The energy of ions in the 
present study is outside this limit and hence, widths calculated 
with Winger’s formula are expected to show worse agreement 
with experiment than some of the other more suitable models 
we considered. 

In an older study, Nikolaev and Dmitriev (ND) ina derived 
a semi-empirical formula for the mean charge states, according 
to which: 


q = Zp[l + iv/Z-v')-^"^]-'^ 


( 11 ) 


where a=0A5, k=0.6, v is the projectile’s velocity, and 
3.6x10^cm/sec. Eq. (11) is designed to reproduce data with 
Zp > 20, at energies of 5 to 200 MeV, mainly on C targets. 
Deviations from experiment on the q values do not exceed 5%, 
according to ca. In the same study, an improved formula for 
the distribution widths optimized for solid ion beam strippers is 
also presented: 


d = do[q[l-iq/Zpy^'^]]^/^ 


( 12 ) 


where do=0.5 and k=0.6. Eq. © is expected to be more accu¬ 
rate for Zp <37 at energies above 20 MeV, as discussed in HI. 
Based on the ND model, E. Baron et al. ifTTl developed their 
own empirical formula to predict the average charge states q. 


which provides a fair agreement with experimental data for 
heavy ions, up to Uranium, in Carbon and Eormvar foils, at 
energies below 80 MeV E|. K. Shima et al. ifTSl developed a 
model for the mean charge states for a large range of ion species 
(Zp >8) in solid targets with 4< <79 at energies below 6 

MeV/u. According to his study: 

q = Zp[l - exp(-1.25X + 0.32X^ - O.llZ^)] 
x[l - 0.0019(Z, - 6) Vz + 0.00001(4 - 6)^Z], (16) 


where the scaled velocity X is given by Eq. 0- The second 
part in Eq. [^is actually a correction term for the non carbon 
solid targets. This model reproduces the experimental data of 
Zp >14 with an agreement of YqjZp <0.04. A reformulation 
of the ND relation is given by the model of To and Drouin 
ca, aiming to reproduce lighter elements (B to Ne) at energies 
up to 7 MeV on C targets. In this model the average charge 
state formula is expressed as: 

q = Zp[l - expi-v/v'z!^p^^)]. (17) 


Even though the formulas in the work of To and Drouin were 
not formulated to reproduce heavy ion/target combinations, 
their work was included in the current study since it is a deriva¬ 
tive of the generally successful Nikolaev and Dmitriev model. 
Einally, A. Leon et al. based on Baron’s work IfTTl . re¬ 
formulated the mean charge state expression by multiplying 
Eq. (13) with a suitable correction factor g’(Z^,Zp). This led 
to improved fits of experimental data at energies of 18 MeV/u 
< Ep < 44 MeV/u for heavy ions with 36< Zp <92 in various 
solid targets (4< Z^ <79).The multiplicative correction factor is 
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given by: 

gf(Zt,Zp) = [(0.929 + 0.269exp(-0.160Z,)) + 

+(0.022 - 0.249exp(-0.322Z,))^] (18) 

where Vp is the projectile’s velocity. 


1.3. Modelling of the equilibrium charge-state distributions 

To succesfully calculate charge state fractions, apart from the 
mean charge state and the width, a model of the shape of the 
distribution is also needed. At low and intermediate projectile 
velocities in light gaseous and solid targets (Z^ <7), the equilib¬ 
rium CSDs tend to be symmetrical. For those cases, the charge 
fractions can be calculated using Gaussian functions: 

Fq = (d V^)-‘ exp[-{q - qfl2d^] (19) 

where d is the distribution width, q is the charge of each state 
(an integer number), and q is the mean charge state of the distri¬ 
bution (in general a real number). For higher projectile veloc¬ 
ities in the same targets, typically for cases in which the mean 
charge state is very close to the Z^, the CSDs become asym¬ 
metrical (these asymmetries are explicitly dependent on the ve¬ 
locity and are not related with the shell effects that we will dis¬ 
cuss next). For those distributions, Baudinet-Robinet et al. ED 
proposed a distribution function extracted from a reduced 
distribution: 

Ft = [2^/^F(y/2)]-if/2-i^-'/2 (20) 


where the chi-squared variable is connected to the charge q, 
mean charge q, and width d as: t=c(Zp-r2-q), c=2(Zp-\-2-q)ld^ 
and y=c(Zp-\-2 - q). 

The CSDs of heavy projectiles, especially in heavy gases 
and solids are asymmetric mainly due to atomic shell struc¬ 
ture. Those cases are far from Gaussian or x'^ distributions. 
K. Shima et al. I22l . proposed the composite of two Gaussian 
functions with the same centroid but different standard devia¬ 
tions (widths) to emulate the asymetrical distributions of exper¬ 
iment. The now different left and right widths are associated 
with the different atomic shells. The usage of such a function 
in the calculation of CSDs requires a model for the estimation 
of double widths which hasn’t been developed so far. Another 
approach for treating the shell effects has been proposed by R. 
O. Sayer 1^ who introduced a modified Gaussian distribution: 

Fq = F„ exp[-0.5t^/(l + et)] (21) 


where t=(q - qo)/p and F^ is the fraction of the most intense 
charge state qo. The shell effects can be satisfactorily repro¬ 
duced (even if they are not explicitly taken into account in Eq. 


21) if the proper values for p and e are chosen. These values can 


be extracted by fitting a suitable function of Zp and projectile 
velocity, ySc, on experimental data. Since there is not a satisfac¬ 
tory amount of data in the heavy ions region (Zp, Zt >20) the 
agreement of this method to experimental CSDs may be some¬ 
what limited. Furthermore, the F^ values must been known in 
advance. Because of these limitations on Sayer’s and Shima’s 



Figure 1: Schematic overview of the FN Tandem accelerator 
and the Accelarator Mass Spectroscopy beamline at the Nu¬ 
clear Science Laboratory of the University of Notre Dame. The 
measurements in this work made use of the MANTIS Browne- 
Buechner Spectrograph, the Multi-purpose Rotational Scatter¬ 
ing Chamber and the FN Tandem. The Wien filter shown in this 
figure upstream of the scattering chamber was not required and 
was not used in the experiment. 


formulas, Eq. ( p^ and Eq. (^) are generally preferred by the 
community for the calculation of CSDs (within their applica¬ 
tion limits) since they can be directly be combined with a large 
variety of semi-empirical formulas for the mean charge state 
q and the distribution width d. In this work we followed the 
same logic and decided to leave out the works of Sayer and 
Shima I 23 II.II 22 I from our comparison to experiment. Hence, 
for the modeling of the charge-state distributions we used only 
the Gaussian and the Baudinet-Robinet reduceddistribution 
shapes described above. 


2. Experimental Procedure and Data Analysis 

The experiment took place at the Nuclear Science Labora¬ 
tory (NSL) of the University of Notre Dame. The incident 
ion beams were accelerated by the 1IMV FN Tandem Van de 
Graaff accelerator. The accelerator mass spectrometry (AMS) 
beamline guided them into the Multi-purpose Rotational Scat¬ 
tering Chamber at the object of the MANTIS spectrograph (Fig. 
[T]). For our measurements the spectrograph was set at a 0° angle 
with respect to the beam axis. Inside the scattering chamber a 
Faraday cup and three 1 pm thick Mo foil targets, were mounted 
on a movable metallic frame with five target positions. The fifth 
position was left blank so that the beam could pass through the 
chamber without interacting with the Mo foil. This setting was 
used during beam tuning. After the target ladder, the beam ions 
could enter the spectrograph where the different charge states 
of the beam could be separated by the magnetic field. 
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Figure 2: Trajectories of the different charge states of the same 
element in the spectrograph magnet. The PGAC detector is po¬ 
sition sensitive along the y-axis. The charge states are sepa¬ 
rated from each other by their charge-dependent position on the 
PGAC detector. The limit in the detector’s active region allows 
the measurement of only those charge states that have a suitable 
trajectory radius in the magnet. 



Figure 3: Charge-states of Cu taken with the ^^Cu beam. The 
spectrum corresponds to a single magnetic field setting. The 
channel number represents position along y-axis in the PGAC 
detector (0 cm to 46 cm). The intensity of each charge state is 
proportional to each peak’s area. 


A Parallel Grid Avalanche Counter (PGAC) detector, with 
active region of 46 x 10 cm, was mounted on a set of rails on 
the top of the spectrograph magnet. At this point the bent beam 
leaves the magnet vertically as it is shown in Fig. By using 
the PGAC it was possible to measure a number of charge states 
simultaneously since the detector is position sensitive. A more 
detailed description of the detector’s operation can be found in 
previous studies of D. Robertson et al. 1241 [25]| . 

The full charge-state distribution for each beam would not fit 
inside the magnetic spectrograph’s acceptance in a single mag¬ 
netic field setting. Therefore, a gradual step-by-step increase 
of the applied magnetic field was used to scan all detectable 
charge states using the detector’s active region. In this process, 
by changing the field from lower to higher values the charge- 
states would be registered by the PGAC detector from higher to 
lower charge, since: 


1^1 cx 4 


magnet 


( 22 ) 


where B is the magnetic field, Imagnet is the current supplied to 
the spectrograph magnet, and a given trajectory radius through 
the magnet is assumed. Its important to mention that no fo¬ 
cusing elements were used before or after the dipole spectro¬ 
graph; so all charge states were transported to the detector with 
the same focusing characteristics as the beam along the disper¬ 
sive direction of the magnet. The detector’s length allowed the 
measurement of only 4 to 6 charge states at each step. As can 
be seen in Fig. |^the charge states appeared in the spectra as 
different peaks along the horizontal axis which represents the 
Y-position inside the PGAC. Scattering in the Mo foil is a con¬ 
tributing factor in the observed widths of the peaks. 

As mentioned before, because of the large number of charge 
states in each distribution, the measurements were performed in 
steps. Between two consecutive steps some charge states were 


chosen to be used as common references allowing the cross¬ 
normalization of the intensities of all charge states along the 
various spectra. To deduce the charge-state distribution from 
the data the following procedure was used: For each charge- 
state q we defined the relative fraction R^, calculated with re¬ 
spect to a reference state (\ref. • 

= A (23) 

qref 

where \q is the intensity of the state q and is the intensity of 
the reference charge state. By using the common charge states 
between adjacent steps as references we were able to calculate 
(for each distribution) all the relative fractions with respect to 
a single state. Having normalized the relative fractions in this 
way, the net fractions F^ could be extracted by: 



(24) 


where the sum is over the total number of charge states in the 
distribution. The uncertainties of the fractions were calculated 
from the statistical errors in peak integration and taking into 
account any overlap of adjacent peaks. 

Having deduced the fractions of all the charge states, the 
mean charge of each distribution was calculated by: 


q 


(25) 


and the distribution widths (assuming Gaussian distributions) 
by: 


q 


( 26 ) 
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For each distribution measured in this experiment the regis¬ 
tered charge states did not constitute 100% of the whole CSD. 
The very low intensity charge states fell below the minimum 
detection limit of the experimental set-up. Moving toward the 
higher charge states the intensity drops very fast since the elec¬ 
tron loss cross section changes significantly; due to the finite 
dispersion of the magnet the high charge states (low m/Q ra¬ 
tio) are expected to be closer to each other, and so, significantly 
more overlapped. These factors limited our ability to measure 
higher charge states in a reasonable time. Regarding the lower 
charge states, no reliable data could be taken for charge states 
with < 10“^% due to the post-foil interactions of the beam 
ions with the residual gas in the beamline vacuum. In such in¬ 
teractions it is expect that the electron capture cross sections 
are higher than the electron loss since the mean charge states 
after the Mo foil will tend to decrease in the residual air. Even 
by assuming -according to our calculations for the experimen¬ 
tal beamline used- that only 0.1% of the beam ions will inter¬ 
act with the residual gas, the effect on the intensity of the low 
charge states with F^ < 10“^ % is still significant due to the elec¬ 
tron captures on higher charge states (especially on those with 
F^ > 5%). These effects are hard to be estimated without fur¬ 
ther measurements or knowledge of the relevant cross sections. 

The significance of the systematic uncertainty induced to the 
charge state fractions due to the missed charge states was es¬ 
timated through a sensitivity test. Three more charge-states 
were added in each CSD and the effects on the charge frac¬ 
tions were calculated. The additional charge states had frac¬ 
tions F^ = 10“^% in order to maximize the effects. The charge 
state fractions varied as a result of this procedure by a factor 
which was found to be at most 0.09%. Statistical errors coming 
from other contributing factors such as peak integration, fluctu¬ 
ated within a range of 1% to 41%. As a result, the error due to 
the missed charge states was considered negligible. 

During the measurements no m/Q interferences from sec¬ 
ondary particles (reaction products) were present. All the nu¬ 
clear reaction channels from the interaction of the beams we 
used (Ni, Co, and Cu) with the Mo were found to have thresh¬ 
olds above 132 MeV while the beam energies we used were 
up to 125 MeV. Furthermore, contaminations in the beam from 
nuclear reactions with the carbon foil in the accelerator’s ter¬ 
minal were eliminated by the 90 degree beam analysis mag¬ 
net. Uncertainties due to possible pile-up were also negligible. 
Throughout the measurements the count rate in the detector was 
monitored and was found to be in no case more than of the order 
of 5000 counts/sec. These rates are too small to induce pile-up 
since the signal processing time in the counter detector used is 
of the order of a few microseconds for each particle. 

3. Results and Discussion 

In Fig. I^the CSDs measured in the current study are pre¬ 
sented. Furthermore, Table [T] includes the fractions of all the 
measured charge states in detail. The statistical error on the 
charge state fractions fluctuated between 1% and 6% in most 
of the charge states. The largest uncertainties were observed 
for the 22-r (~ 17%), 23+ (~ 30%), and 24+ (~ 41%) charge 


states of Ni, Co, and Cu respectively. The increased errors at 
the highest charge states are due to the low statistics of the cor¬ 
responding peaks in combination with the overlap of adjacent 
peaks for these higher charge states. In order to be consistent 
with literature data, we assigned the CSDs to the emerging en¬ 
ergies -i.e. for ions that have already lost energy by traversing 
the thickness of the target. In this way there is a direct cor¬ 
respondence between equilibrium CSD and projectile velocity 
(which is not the case when the initial beam energy is assigned). 
The energy loss was calculated using SRIM-2013 1^ and was 
found to be approximately 20 MeV for all three beams. The 
errors introduced to the mean charge state due to the energy 
loss calculation with SRIM were estimated to be less than 0.01 
charge units. 

All three distributions were found to be asymmetric due to 
shell effects. This asymmetry is evident by examining the ra¬ 
tios of the fractions for each CSD. In the case of a symmet¬ 
ric (Gaussian-like) equilibrium CSD, the logarithm of the ra¬ 
tios Fg+i/Fg is linearly varying with the charge q. This is be¬ 
cause, at equilibrium, the ratios cr^,^+i/cr^+i,^ (where cr^,^+i is 
the electron loss cross-section at the charge-state q and cr^+i,^ 
is the electron capture cross-section at q-rl) are approximately 
proportional to e“^ and Fq+i/Fq=crg,q+i/crg+i g ||Z7]| . The above 
statement implies that the single electron exchanges are domi¬ 
nant. However, when shell transitions or multiple-electron ex¬ 
changes occur the relationship between ln(F^+i/F^) and q is not 
linear, resulting in asymmetrical CSDs. In Fig. the discon¬ 
tinuities at q=16-h, q=17-r, and q=18-r (in Co, Ni, and Cu re¬ 
spectively) are consistent with the L-M shell transitions as can 
be shown by examining the corresponding electronic configu¬ 
rations ||28l. The fractions of charge states that correspond to 
closed shells, or sub-shells, are significantly enhanced resulting 
in ’’kinks” on the plot. For the case of Co, the discontinuity at 
q=22+ seems to correspond to the 2s-2p sub-shell transition. 

In Table all the calculations of the mean charge q using 
various semi-empirical models are presented. The experimen¬ 
tal values have been calculated by Eq. (4). It is important to 
note that the mean charge of the distribution doesn’t have to 
be an integer number. The maximum charge state (the most 
intense) is the integer which is closer to the mean value. In 
terms of their predictive power for the maximum charge state, 
Schiwietz et al. and Shima et al. offer the best agreement to 
experimental data. The agreement of the other models exam¬ 
ined in this work is within ±2 charge units. The model by To 
and Drouin provided a better than expected agreement with ex¬ 
periment considering that it is designed to reproduce CSDs of 
lighter elements. On the other hand, the formula of Leon et 
al. (devised for heavier elements and higher energies) had the 
worst agreement of all models to our measurements. Regarding 
the mean value. Winger’s model is in agreement to the experi¬ 
mental values within the uncertainty limits in most cases. 

In Table the predictions of the semi-empirical models for 
the distribution widths are presented. The experimental widths 
have been calculated using Eq.(5). To compare each width for¬ 
mulation’s performance to the others, the experimentally deter¬ 
mined mean charge states q were used in all the calculations. 
By neglecting shell effects, most of the models diverged from 
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Figure 4: Charge state distributions of ^^Ni (black inverted tri¬ 
angles and line), ^^Co (red circles and line), and ^^Cu (blue up¬ 
right triangle and line) beams in a lyum thick Mo foil. No error 
bars are visible for statistical errors less than 3% due to the size 
of the point markers. The experimentally deduced mean charge 
states for each distribution are ^=18.45 for Ni, ^=18.20 for Co, 
and ^=19.33 for Cu. The most intense charge states were 18-r, 
18-r, and 19+ for Ni, Co, and Cu respectively. 



Figure 5: Logarithm of the ratio of two adjacent charge state 
fractions as function of the charge of ions detected with the 
MANTIS spectrograph in the present study. Lines are used to 
guide the eye. Black inverted traingles and a black line corre¬ 
spond to Ni ions on Mo using the ^^Ni beam, red circles and line 
correspond to Co ions on Mo using the ^^Co beam, and blue up¬ 
right triangles and line correspond to Cu ions on Mo using the 
^^Co beam. Divergence from a straight line is attributed to shell 
effects and multiple-electron exchange processes. All observed 
distributions demonstrated such deviations which correspond to 
non-symmetrical CSDs. 



Figure 6: Equilibrium CSD of Cu in Mo, as measured in the 
present study (blue line, upright triangles) and as reported by K 
Shima et al 1291 (magenta line inverted triangles). The distribu¬ 
tions are assigned to emerging projectile energies (after energy 
loss in the Mo foil). Mostly excellent agreement between the 2 
measurements is observed. 


the experimental widths. A systematic overestimation was ob¬ 
served on the values calculated with Nikolaev-Dmitriev (ND) 
and Betz models, while the Baron et al. and Winger et al. mod¬ 
els systematically underestimated the distribution widths. The 
model by Schiwietz et al. was shown to be in better agree¬ 
ment with experiment in reproducing the widths of the mea¬ 
sured CSDs with very small deviations. This success we at¬ 
tribute to the fact that this model was based on fitted data that 
included a significant number of asymmetrical CSDs. 

The only available data in literature, relevant to the current 
work, are coming from a previous study by K. Shima et al. 1291 . 
Shima et al. used a ^^Cu beam impinging on a Mo foil. The re¬ 
sults are plotted in Fig. [^in comparison to data from the present 
work. Comparing the two data sets we see a mostly excellent 
agreement to each other considering the small energy difference 
of the two cases. The fractions of the charge states 19+ and 20+ 
are identical within the statistical errors while deviations are ob¬ 
served on charge-states with fractions lower than 15%. In terms 
of the mean charge state and the distribution width (see Tables 
[^and|^, the deviations are within the statistical errors. 

3.1. A combinatoriaV’ prescription for reproduction of the 
charge state distributions 

The results of this study support a systematically better 
agreement with experimental data of the formulation of Schi¬ 
wietz et al. oa for the distribution widths and of Winger et 
al. ca for the mean charge states q. Combining these models 
with a Gaussian or a reduced distribution function, we ex¬ 
plored a more realistic way to reproduce the shape of CSDs in 
the region of our experimental study (Z^ ~28 and Zp=42). We 
compared this “combinatorial” model to calculations utilizing 
a single model to estimate both the distribution width d and the 
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mean charge state q. Each of these calculations was performed 
in two variants; one assuming a Gaussian shaped charge state 
distribution and one assuming a Baudinet-Robinet type reduced 
distribution. The results are presented in Fig. [^(a-f) while 
on Table 1^ the values of the calculated CSDs are presented 

as extracted from chi-square goodness of fit tests. 

From the comparison in Fig. |7]it is suggested that the com¬ 
binatorial model combined with a reduced distribution func¬ 
tion produces a qualitatively better description of the experi¬ 
mental data especially in the regions of low intensity charge 
states near and at the tails of each distribution. At the higher in¬ 
tensity charge states around the mean no significant difference 
between the reduce and Gaussian-shape charge distributions 
is observed. Nevertheless, the overall agreement offered by the 
combinatorial model (combined with any of the two distribu¬ 
tion functions we discussed) is improved in comparison to the 
other models as demonstrated by the lower on average value of 
X^ (see last row of table |^. 


4. Conclusion 

Equilibrium charge state distributions of ^^Co, ^^Ni, and 
^^Cu beams passing through a lyum thick Mo foil have been 
measured. A variety of semi-empirical models for the mean 
charge state q and the distribution width d of equilibrium charge 
state distributions were compared with our experimentally de¬ 
termined charge state distributions. Furthermore, our study sug¬ 
gests that an improved agreement of the calculated equilibrium 
CSDs to the experimental data in the region of study (Z^ ~28, 
Zp=42, and 2 MeV/u) can be obtained by using a com¬ 
bination of models to describe the equilibrium CSDs. In this 
’’combinatorial” description, the formulation of Winger et al. 
ca is used to calculate the mean charge state q, the work of 
Schiwietz et al. oa is followed in the calculation of the dis¬ 
tribution width J, and a reduced x'^ or a Gaussian function can 
be used to describe the shape of the charge state distribution. 
Despite the improved agreement of this ’’combinatorial” model 
it is still a phenomenological prescription with obvious limita¬ 
tions. Realistic first-principle based simulations of equilibrium 
and non-equilibrium CSDs of heavy ions would be the ideal 
way to go forward. Until this is reliably possible, it would be 
beneficial to extend the current set of experimental data to cover 
the 19< Zp, Zt <54 region at various energies below 10 MeV/u. 
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(C) 



(d) 




(e) (f) 

Figure 7: Calculated CSDs using different semi-empirical models including our “combinatorial” formulation in comparison with 
experimental data. The left column of the figure shows calculations assuming a reduced type shape of the distribution, while 
the right column shows calculations using a Gaussian shape. Figures (a,b): Comparison of experimental data from this work to 
calculations for a ^^Co beam. Figures(c,d): the same for ^^Ni. Figures (e,f): The same for ^^Cu. In all figures, experimental data 
are represented by black stars and line, the calculations using the work of Schiwietz et al (S in figure) are presented with red circles 
and line, the calculations using the work of Winger et al. (W in figure) are presented by magenta inverted triangles and line, and 
the combinatorial formulation introduced in this work (W-rS) is presented with blue upright triangles and line. Values of;^^ for the 
comparison are presented in table 9 





































































































































































Table 1: Fractions percent (%) of all the measured charge states in the present study. Each row of the table corresponds to a 
different beam in Mo. 


Ion 

E(MeV/u) 

Fi4 

Fi5 

Fi6 

Fi7 

Fi8 

Fi9 

F20 

F21 

F22 

F23 

F24 


1.769±0.005 

0.38±0.02 

1.67 ±0.08 

5.71±0.30 

21.55±0.37 

30.24±0.42 

24.98±0.22 

11.86±0.14 

3.04±0.08 

0.49±0.02 

0.08±0.02 

- 

60Ni 

1.510 ±0.009 

0.35±0.02 

1.58±0.05 

4.82±0.12 

13.97 ±0.33 

31.04±0.25 

28.35±0.51 

14.18±0.21 

4.71±0.16 

0.99±0.17 

- 

- 

63Cu 

1.783 ±0.018 

- 

0.52±0.02 

1.99±0.05 

5.96±0.07 

14.78±0.21 

31.74±0.71 

27.06±0.64 

13.23±0.23 

3.88±0.09 

0.74±0.02 

0.10±0.04 


Table 2: Comparison of experimentally and theoretically determined values of the mean charge state q in Mo. Predictions of the 
various semi-empirical models considered in this work are shown. Each row of the table corresponds to a different beam. The last 
row gives for reference the experimental value from literature for Cu which is in agreement with our own measurement. 


Ion 

E (MeV/u) 

Exper. 

Schiwietz [13J 

Shima 1221 

Baron [VT\ 

Nik.-Dmit. 11161 

Drouin [19] 

Winger Il51 

Leon 1201 


1.769 +0.005 

18.21 +0.12 

18.02 

17.73 

17.67 

19.12 

18.55 

18.37 

16.59 

60Ni 

1.510+0.009 

18.45 +0.14 

17.90 

17.57 

17.32 

18.96 

18.25 

18.49 

16.26 

®Cu 

1.783 +0.018 

19.34 +0.20 

19.12 

18.80 

18.70 

20.31 

19.66 

19.57 

17.52 

®Cu 

1 A . 

1.767^ 

19.3E 









^ As reported by K. Shima et al. HU 


Table 3: Comparison of experimentally and theoretically determined values of the distribution width d in Mo. Predictions of the 
various semi-empirical models considered in this work are shown. Each row of the table corresponds to a different beam. The last 
row gives for reference the experimental value from literature for Cu which is in agreement with our own measurement. In order to 
properly compare the predictions of the models with the experimental distribution width, the experimental value of q was used as 
common input in all calculations. 


Ion 

Target 

Experimental 

Schiwietz 1141 

Baron ri71 

Nik.-Dmit. [1C| 

Winger [15| 

Betz (4] 


Mo (foil) 

1.31 +0.02 

1.34 

1.23 

1.48 

1.19 

1.40 

“Ni 

Mo (foil) 

1.33 +0.02 

1.37 

1.26 

1.52 

1.20 

1.43 

«Cu 

“Cu 

Mo (foil) 
Mo (foil)i 

1.35 +0.02 
I. 33 I 

1.38 

1.28 

1.54 

1.22 

1.45 


^ As reported by K. Shima et al. 1291 


Table 4: values of the calculated CSDs. In the calculations were used either combination of models or single models for the q 

and d. For explanation of labels see text and figure(W): Winger et al., (W+S): this work, (S): Schiwietz et al., (G): Gaussian 
distribution, {x^)\ Reduceddistribution. 



W-hS (G) 

W(G) 

S(G) 

w+s Or") 

wor") 

sor") 


2.48 

7.22 

2.58 

3.85 

11.31 

6.17 

“Ni 

1.74 

7.20 

17.85 

1.81 

5.92 

22.08 

“Cu 

5.33 

16.45 

4.12 

4.66 

9.85 

4.93 

Average 

3.18 

10.29 

8.18 

3.44 

9.02 

11.06 
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